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ABSTRACT 


This thesis describes the conversion of a computer 
program from Fortran IV used by the CDC computer to Fortran 
IV compatibie with the Naval Postgraduate School IBM 3033 
systen. The converted program, called TURB2, estimates the 
Magnitude of the refractive index structure parameter, 
Сп2 (2), for a dry atmosphere in horizontal layers. The 
altitudes of the layers depend upon the corresponding alti- 
tudes of conventional meteorological rawinsonde bailoon 
data. The data input is a formatted file called  TURB2 
DATAIN and the output consists of the value of dry Cn? at 


each altitude along with a graph of log Cn? vs. altitude. 
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I. INTRODUCTION 

Many areas of modern technology are concerned with the 
propagation and distortion of  eiectro-magnetic radiation in 
the atmosphere. In particular, atmospheric optical turbu- 
lence seriously degrades visible light astronony, satellite 
photography and laser propagation. 

Studying and attempting to predict atmospheric optical 
index of refraction fluctuations has been going on fora 
number of years [Ref. 1]. With the ascendancy of modern 
computer methods during the last twenty years, the 
predictions have grown significantly in accuracy and commen- 
surately in complexity. Н.Е. Hufnagel published a model 
[ Ref. 2], in 1974 based on empirical data collected at 
astrononical observatories. However, it did not consider 
atnospheric dynamics other than an average wind speed. In 
1577, T.E.  VanZandt published a nev model (Ref. 3]. This 
model is the most comprehensive and accurate one developed 
to date and attempts to include actual temperature and wind 
shear gradients that force the production of optical 


turbulence. 


A. BACKGROUND 


The computer code based on  VanZandt's model was written 
by J. Warnock at the National Oceanographic and Atmospheric 
Administration's Atmospheric Dynamics Aeronomy Laboratory in 
Boulder, Co. The latest version of the program, that is 
constantly being updated, has a modular format and consists 
of twenty subroutines and functions. The code was written 
to be executed on NOAA's CDC Cyber 750 computer using 


Fortran Iv: The primary concern of this thesis was the 


translation of the subroutines as received to code that 
would be compatible with the Fortran IV used at the Naval 
postgraduate School's (NPS) ІВМ 3033 computer system and to 
exercise this program on specific atmospheric wind and 


temperature profiles. Additionally: 


1. a main program was written to consolidate the subrou- 
tines, read the input data file and print the tabular 


output, 


2. a data processing subroutine was written to transform 
the raw input data into the variables required in the 


following subroutines, and 


3. a Disspla routine was written to produce a graphical 
presentation of the computer model output and to permit 


comparison with whatever reference values were input. 


B. GOALS 


With the recent interest that the Department of Defense 
has shown in visible wavelength laser systems, an accurate 
computer model for atmospheric optical turbulence is essen- 
tial for meaningful systems analysis. Idealiy, this model 
should be relatively fast and easy to use and should be able 
to operate with readily available data. 

The program that has been developed is capable of satis- 
fying these requirements. The raw input data comes fron 
conventional rawinsonde printouts. These are available in 
the computer record archives of the National Weather Service 
for hundreds of stations, ashore and afloat, and going Lack 
quite a number of years. The program produces tabular 
and/or graphical output from the rawinsonde data entered in 
the data input file (hereafter DATAIN file). For a given 
location and date, the program estimates the turbulence 


Structure parameter, Cn?, between each input data altitude. 
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Because it is rather easy to use, it is hoped that the 
program will prove to be beneficial for physical meteorology 
and laser technology classes. It is also very instructive 
to modify the input data to determine the effects of several 
atmospheric variables upon the turbulence structure paran- 
eter. For example, it is possible to determine the effect 
of variations in the temperature or wind shear profiles upon 
Ene final output. 

The program should be useful as a reference to 
researchers who are trying to measure the Cn2 profile exper- 
imentally. Assuming the availability of a nearby rawinsonde 
launch site (within 50 to 100 miles), the computer model 
will produce Cn? estimates for comparison with experimental 
data. Ideally the experimental measurement device or the 
model will be fine-tuned as a result of detailed comparative 
analyses. 

Finally, the program also has the potential to be usei 
aS a subroutine to estimate the Cn? profile for other 
programs. For example, the G.U.7.S. program presently used 
in AE-4706 makes use of a rather crude approximaticn for 
Gne. Using the much more realistic values obtained fron 
this computer model would serve to improve tae accuracy of 
ες G.U.T.5. output. 

The primary goal of my work was to translate the CDC 


code into a smoothly running and error-free program for the 


IBM system at NPS. Also, the format of the program as 
received, was somewhat rough and required a significant 
amount of time to understand and be able to use. AS a 


secondary goal then, it was desired to make the program as 
"user friendly' as possible so that it could be used for 
class projects and by researchers with a minimum of instruc- 
twon and difficulty. Finally, since the Cn? parameter can 
vary erratically with altitude, it was thought that sinple 
tabulated output would be of limited utility, especially as 


11 


an instructional aid. Therefore a graphical 
capability was desired. Additionally, it is 
that a researcher using the program would find 
be able to compare several profiles, perhaps 


different methods. Thus the capability to read 


presentation 
quite likely 
it useful to 

obtained by 


Several data 


files and plot the ccmparative data was also incorporated in 


the program. 
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II. APPROACH TO PROBLEM 

The model program was received from NOAA in printout 
form and was typed onto my 'A' disk on the NPS system with 
each subroutine filed as a separate Fortran program. During 
this phase, it became evident that the full 60 bit capa- 
bility of the CDC computer had been utilized . To achieve a 
comparable level of accuracy on the IBM system, it was 
necessary to use double precision. Rather than change all 
the variables to double precision before more familiarity 
was gained with the program, it was decided to use the 
autodbl function which would automatically change all the 
variables and functions to IBM doubie precision. MI 
Hilleary at the computer center was most helpful in this 
endeavor. 

Because of the complexity of the program and its modular 
structure, it was decided that a ground-up module test 
method would prove most advantageous. Thus once all the 
Subroutines were copied and checked, a short WATFIV progran 
WaS written around each subroutine to test its correct oper- 
асот. The WATFIV compiler was utilized because of its 
excellent error messages. mac program supplied the 
required input variables, called the respective subroutine 
and displayed the desired variables which were spot-checked 
using a handheld calculator. During the initial parts of 
this phase, Lo attempt was made to integrate the subrou- 
tines, i.e. each program tested one subroutine or function 
only, in an attempt to minimize the amount of debugging 
required. However, once all the subroutines had been tested 
individually, some integrated testing was conducted. The 
modular testing procedure additionally served as an excel- 


lent means of gaining familiarity with each part of the 


ШЕ 


program. This proved to be of great benefit in later stages 
of debugging and modification. 

Simultaneous with the testing and debugging process, it 
was realized that actual rawinsonde data would be required 
once the program was ready to run. To avoid a lag in the 
development, an effort was made to procure some representa- 
tive data fairly early. The dates and locations were chosen 
to coincide with a detailed experiment run by the University 
of Arizona and the Air Force [ Ref. 4] at Sunset Canyon, AZ. 
and Truckee, Neve. іп 1977-78. The Naval Environmental 
Prediction and Research Facility (NEPRF), located іп 
Monterey, Ca., provided rawinsonde data upon letter request. 
The rawinsonde data output was in the form of raw tabulated 
values for the mandatory and significant reporting levels 
and a summary table which listed the desired variables 
including atmospheric pressure, temperature, humidity, wind 
speed and direction and absolute altitude. Additionally 
graphs were also received but were not utilized. Тһе sumna- 
rized values were filed in the proper format in the DATAIN 
file. 

The next step was to write a data processing subroutine 
that would: 1. process the raw data into the units required 
by the subroutines and, 2. calculate the necessary derived 
variables such as the humidity gradient and the wind shear. 

With the completion of the data processing subroutine, 
all was in readiness for the development of the main 
program. The functicn of the main program was to step the 
calculations from the first data input set (corresponding to 
the lowest slab of a given atmospnere) sequentially througn 
each slab until the last data set was reached. The main 
program also initialized a number of common variables that 
were used in various subroutines throughout the program. 
Lastly, it served as a framework for the proper sequencing 


of the subroutine calls. 
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As night be expected in complex programs, even with the 
extensive modular testing previously conducted, a signifi- 
cant amount of debugging was required once the program was 
finally integrated. When attempts to decipher the Fortran 
error messages proved fruitless, the WATFIV compiler was 
again utilized. The main difficulty in debugging was trying 
to trace the errors through a series of  imbedded subrou- 
tines. For example one of the most difficult problems to 
detect at this point was caused by an underflow error. Ihe 
EXPNX function was written to prevent underflow or overflow 
errors from occurring. As received, it tested to see if 
variables exceeded the limits of the CDC computer (approxi- 
mateiy 1.£+300). During the early testing, it was realized 
that these limits would have to be changed to correspond to 
IBM system 303315 capacity (approximately 1.5+75). If the 
variable in question exceeded these limits, the function 
EXPNX would revalue the variable to the respective limit 
imposed by the computer. What was not immediately apparent 
however, was that the corrected variable was subsequently 
multiplied by another small number (on the order of 1. 5-3) 
in another subroutine which uitimately caused the underflow. 
While in itself, the error was relatively easy to detect, it 
proved difficuit and time consuming to trace since the 
offending subroutine was buried in other subroutines. 

Che technique that also proved effective was to drawa 
large flowchart which listed each subroutine and function 
contained in the progran. Additionally eacn subroutine or 
function was annotated with the other subroutires that 
called it as weil as those that were called by it. This was 
ΙΓ; 1. Gaining the Thig picture! of the struc- 
ture of the program, 2 Checking that the reguired vari- 
ables were indeed entered or printed and that each was 
listed in the call statement, 3. detecting several yari- 


abies that had inadvertently undergone name changes in the 


ПР 


genesis of the program at NOAA and had subsequently become 
sources of confusion in the translation, and 4. tracing 


suspect variables back to their origin. 
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III. SOLUTION 


А. THEORY 


VanZandt et al, [Ref. 3], assume the existence of thin 
Lorizontal layers of turbulence in an otherwise  non- 
turbulent atmosphere. These layers are evident in radar 
[Ref. 3] and direct balloon measured Cn? profiles, (Ref. 5] 
and [Ref. 6]. They vary with time and location, though they 
are generally of limited horizontal area. The large scaie 
turbulent eddies in this thin region are inherently inhomo- 
geneous, anisotropic and non-steady. Nevertheless, within 
the turbulent layers, the turbulence approaches homogeneous 
and isotropic conditions if the scale lengtn is taken small 
enough (i.e. scale lengths much shorter than the thickness 
of the layer). Furthermore, it is assumed that all the 
different sized eddies come into existence within a short 
time span, relative to the life span of the turbulent eddy 
itself. With these approximations, the velocity fluctua- 
tions can be considered to be steady-state and the turbu- 
lence structure parameter, Cn?, can be estimated for the 
turbulent layer. To estimate the turbulence throughout an 
entire vertical shaft of atmosphere, it is also necessary to 
estimate the fraction, F, of the entire shaft tnat is turbu- 
lent. The calculation and combination of these two esti- 
mates form the basis for this model and the program. 

ror the case o£ homogeneous, isotropic turbulence, 


Tatarski [Ref. 1], estimates Cn? by; 


Cn2 = ἂ2 ο αἲ ο 1, 4/3 e M? (3.1) 
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where a is a universal constant generally taken to be 2.8, 
о" 15 the ratio between the eddy viscosity and the eddy 
thermal diffusivity generally taken to be 1.0, Le is the 
outer scale length of the turbulence and M is the vertical 
gradient of the refractive index. 

Le is assumed to be of the order of the thickness of the 
turbulent layer itself. An expression for le, also given by 


Tatarski, is; 
Lem (e 7 Sie (3299 


where e is the dissipation rate of turbulent kinetic energy 
and S 15 the vertical gradient in the horizontal wind 
velocity (i.e. DEA 2 Z) - Rather than try to measure these 
values, the developers of this model selected a value of Le 
that would normalize the theoretical Cn? to that produced by 
their experimental radar measurenents. This value of Lə 
turned out to be 10 meters. 


The expression for M is given by [Β6Ε. 31 85; 


M = —77.6 e 1.E-5 e ( p e 91n@8 / Т92) (3239 
° ( 1 + 15500 g/T e[ 1-.5( 21n gZ22) Z (Ə1n@ / 32) ] 


where p is the atmospheric pressure in millibars, T is the 
absolute temperature in Kelvin, O is the potential tempera- 
ture in Kelvin, gq is the specific humidity and z is the 
altitude in meters. To calculate М, the vertical resolution 
of the meteorological data sSnould be on the order of the 
layer thickness, Lo. Current rawinsonde measurements are 
somewhat coarser, on the order of 100 meters (for routine 
equipment, special equipaent is available which is capable 
οἱ producing a much finer dataset), however the rawinsonde 
data is freguently the best data available. Using eqns. 
(3.3) and (3.1), it is possiLle to estimate the value of Cn? 


in the turbulent layers from routine rawinsonde data. 
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A slab of atmosphere on tke order of 10 meters thick 
will usually contain at least one set of turbulent and non- 
turbulent layers. Having previously estimated the value of 


Cn2 in the turbulent layer, the average Cn? for slab i is; 
Cn2 (i) = Cn2(i) ° F(i) (224) 


where F(i) is the mean fraction of slab(i) that is turbu- 
lent. There are no expressions to solve for F directly so a 
statistical distribution is used. 

First one calculates the wind shear, which is the 
difference between the horizontal wind velocity vector at 
two consecutive altitudes (note that any vertical wind 
component, such as due to convective flow, is neglected). 
At any point, the total shear, ‘St, can be thought of as 
being composed of two parts, the mean Shear, Sm, and a fiuc- 


ышат лч pact, Sf. Then at any point; 
St = Sm + S£ (3.5) 


The coarseness of the rawinsonde measurements will provide 
σου accuracy for the calculation of Sm, but will 
preclude the direct calculation of Sf which would resuire 
the separation of data points to be on the order of several 
meters. Thus it was assumed that Sf has a Gaussian distri- 
bution 7( ) and a standard deviation, sigma, which is inde- 


pendent of Sm so that; 
Z (Sf/sigma) - exp( -Sf2/2esigma2) / (2¢pi) 172 (3.6) 


Then the probability, Ё, of a total shear in the increment 
ο τοσο τας is; 


f£(St) dSt - Z[ (St-Sm)/sigma] d(St/sigma) (2-7) 


It was also assumed that a region is turbulent only if St 


exceeds a critical shear, Sc. Then F is given by; 
r [Sel τ 


F(|taucrit],]taumean]) - fc suas: + £ (St)dSt (35,28 
2 сө “|51 
- P(-|]taucrit]-|taumean]|) * Q(l|taucrit|-]taumean|) 
where jtaucrit] = {Scj|/sigma, ltaumean| = |Sm|/sigma, P( ) 
is the cumulative distribution function of Z() and Q( ) = 1 


- P( ). Thus, F varies with |taumean|], the magnitude of the 
wind shear and inversely with |taucrit], the atmospheric 
stability. 

From есп. (3.8), F is a function of two  variarles, 
taucrit and taumean, which in turn are functions of tu 
variables, Sc, Sm and sigma. Sn, as mentioned, is computed 


from the rawinsonde data. Sc is caiculated from; 


5с -( (4 “діІп0/92) / ἃς 11/2 (325% 
- 2(g *31n0/3z 1/2 


Ihe critical value for Richardson's number, Rc, for the 
onset of turbulence is taken to be .25 and this value was 
substituted in the second equality of eqn. (3.9). The value 
of sigma is taken as .01/s in the troposphere and .015/s in 
the stratosphere, [{Ref. 3]. Fith these variables in hand, 
Ғ(і) for each slab can be calculated from есп. (3.8) and 


uitimately Cn? from «ап. (3.4). 


B. METHOD 


The program as received was similar to IBM Fortran IV 
overall, although there were some significant differences. 
The CDC to IBM Conversion Guide, fRef. 7}, proved most 
helpful as did the staff at the Computer Center and local 


expertise in the Aeronautical Engineering Dept. at NPS. the 


2:9 


nost frequent modifications that were made have been listed 
io Takle 1. 


TABLE 1 
CDC to IBM Translation 


Item CDC code IBM translation 
Comment card uses * or C SubSritutrerciror” x 
їп COL, #1 
Continuation uses '1' in substitute * for 1 
card ccl. #6 
Statement uses А-1%В-2%С-3 write each statement 
separation for example on IUE lines, 
16. = 
B=2, etc. 
initialize uses A=B=C=0 write each statement 
Statements ог "η 
different lines, 
l.e. A=0 
B=0, etc. 
Program uses Program delete Program 
statements Statement at the Statement 


start of the main 
subroutine 


Name length permits up to 7 rename variables so 
of variable Characters inh that maximum length 
name doesn't exceed 6 
Characters 
паше теп іа. same as above same as above 
of subroutine 
Delimiters uses PRINIUABC! substitute 
for example τος place, ot T 
Print uses PRINT*A,B,C substitute 
statements or PEINT 100,X Wamoie (Ss, 100) x 
as required 
Integer uses LAST=NQUAD/2 rewrite using 
division where LAST and different approach 


NQUAD are ведена 
and NQUAD is odd 


Once the program was running satisfactorily, the Dissla 
graphics routine was developed. Since the program was 


written as single precision and then converted into double 


precision through the use of the autodbl function, all of 
the output values were necessarily of double precision type 
and it was impossible to selectively override the function 
SO as to declare some of the variables as single precision. 
Because Disspla will not accept double precision, two alter- 
natives presented themselves at this point, either; 1: 
write a separate, single precision plotting  progran that 
would read the main program's data output file (hercafter 
the DATAOUT file), as single precision, and plot the values 
or, 2. convert the entire program to double precision except 
for two single precision variables which would serve as the 
x and y coordinates for the Disspla routine. Since tne 
program for alternative #1 was already working, it seemed 
the most expedient choice and therefore a separate, short 
plotting program was developed. Once it was runLing 
however, it became apparent that having to run first one 
program and then another was overly cumbersone, especially 
for a program to be made available for general use. 
Alternative #2 was then developed, i.e. all the program 
variables were declared as REAL*8 except for two single 
precision, dummy variables, PLOTAX and PLOTAY, and the plot- 
ting routine was incorporated as a section of the main 
program. This proved very effective because the user how 
received the graphics at a Tektronics 618 terminal followed 
by the tabulated data ina DATAOUT file which is available 


for browsing or printing from the user's Filelist (РІЕ5 


C. INITIAL CALCULATIONS 


The main program first declares all the variables and 
initializes the necessary constants which are used in the 
program and reads the program options (i.e. whether or not a 
plot is desired as well as the number of runs desired). 


Next the program reads the number of data sets, N (i.e. the 
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number of levels reported by the rawinsonde), from the 
DATAIN file. The variable format is I5. Following that a 
"Го loop", iterated N times, reads the data sets from the 
DATAIN file. For each level, the data set format is 
eee 10.5), F10.2 and the variables are inthe following 
order: pressure (millibars), temperature (Celsius),  rela- 
tive humidity (%/100), wind speed (meters/ second), wind 
direction (compass heading in degrees relative to true 
North) and absolute altitude (meters). The respective vari- 
abie names are: PRESS, TEMP, HUMID, WS, WD, and Z. 

The data processing subroutine  (DAT2RO) is called with 


these six variables plus N. The subroutine calculates; 
TEMPA (the absolute temperature, TEMP + 273.) 
TZKARY (the thickness of the slab, Z(I*1)-Z(I)) 


DODZRY (the average gradient of the humidity through the 
Slab) 


PARY (the weighted average of the pressure for the slab) 


TARY (the weighted average of the temperature for the 
slab) 


CARY (the weighted average of the humidity for the siab) 


PARY, TARY and QARY were originally taken as simple aver- 
ages, but were changed to weighted averages in order to 
smooth the roughness introduced when point rawinsonde values 


are used for continuously varying (spatially and temporally) 


physical parameters. The generic weighting 2cuation used 
was; 
een) τπτ (hal), + 3. Y (I) (3. 10) 


ο ο νο гы T Y (1+2) ] Z 8. 
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The value of THETA, the potential temperature, is calcu- 


lated next fron; 


THETA (I) = [ (PRESS (1) /PRESS (I)) **(R/Cp) JeTEMPA (3.11) 


where R is the gas constant and Cp is the specific heat. 

The wind vector is then broken down into the usual u and 
v components (UX and UY respectively) and the wind shear is 
calculated from; 


SHEAR(I) = [ (UX(I*1) - UX(I) )2 + (3982) 


e 


( UY(I*1) - UY(I1) )2 1142 / ΤΗΚΛΕΥ͂ (Τ) 


Next, the stability parameter, STBARY, is calculated by; 


STBARY(I) = Gef LOG(THETA (I+1))-LOG(THETA(I)) 1 (3.13) 
ZAR ZEE Zo CT) 


iinally Richardson's number, RIARY, for the slab is 


calculated by; 


ЗТАВҮ(І) - 5ТВАВҮ(І) / (5НЕАЕҮ(ЦТ)) 2 (3. 14) 


The values of TARY, PARY, QARY, STBARY, SHEARY,  DODZRY, 
RIARY and THKARY are ali arrays vitni N- le ΠΤ These 
arrays are returned to the main program. The primary 'Do 
loop" of the main program, which iterates over the altitudes 


corresponding to the input data, then calls the sulroutines 


SERSG M; STBSGM ang 26555 INTLGS calls each of the 
remaining subroutines and functions. A flowchart, Figure 
3.1, depicts the general flow within the progran. Rather 


than describe each subroutine, a brief summary of each is 
included as Appendix B. 
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The end product is the value of PCNSQ (ог Cn?) calcu- 
lated for each slab. The average altitude of the slab, 
MIDALT, is used for the graph and tabulation. The Cn? 
values with the corresponding MIDALT's are stored in two N-1 
element, single precision arrays,  PLOTAX and PLOTAY respec- 
tively. At this point, if the number of runs is specified 
as '2* in the DATAIN file, the program reads in all the 
datasets again and recalculates the Cn? profile (this option 
Will be further discussed in Chap. 4). Next, if specified 
(i.e. if PLOT = "2% іп the DATAIN file), the Disspla routine 
is invoked which plots the value of log Cn? on the x-axis 
vs. MIDALT on the y-axis. The graph displayed is available 
in hardcopy form at the Tektronics 618 terminal. Finally, 
the Cn? and MIDALT values are also listed in the  TURB2 
DATAOUT file which has been transferred to the user's 'А* 


disk. This file may then be browsed or printed as desired. 


D. INPUT FORSAT 


The following section describes the step-by-step proce- 
dure to be followed when executing TURB2. Prior to execu- 
tion of course, it is necessary to obtain copies of the 
reguired programs, namely TURB DATAIN,  TURB2 FORTRAN and 
optionally; REPPLOSOGDATATNS Additionally, an exec, TURB2 
EXEC, is available and iS very convenient to use. 

The first step is to type the rawinsonde data into the 
TURB DATAIN file. The easiest method is to copy, rename and 
write-over an old DATAIN file (see Appendix A), however the 
format is listed in Table 2 . 

When the rawinscnde data is received, each data set 
(i.e. a set of all tne measurable variables at a given alti- 


tude) will fall into one of three categories, either; 


15 the dataset has balloon measured data (pressure, 


temperature and humidity along with radar altitude) tut 
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TABLE 2 
TURB DATAIN File Format 


Card (s) Cols. Type Variable kemarks 
1 20 I5 PLOT enter '1' if a plot is 
desrresd;wt2v"wr not 
2 20 15 #RUNS enter the number of 
runs desired (1 or 2) 
3 162/20 πο Ν number of datasets to 
be entered 
4->N+4 I Ετος” PRESS atmospheric pressure 
(milli bars) 
п 210120 5 TEMP temperature (degrees 
Celsius, not absolute) 
п 215320 EOTS HUMID relative humidity (7) 
£ 31-40 ΙΟ» WS wind speed (meters/sec) 
п 41-50 TE WD wind direction (compass 
direction, degrees) 
и 5 12619 ЕШ 2 absolute altitude 
(neters) 


Note: If #RUNS = 1, the DATAIN file is Complete. If RUNS = 
copy cards 3 thru 0+ 0 and insert them immediately following 
card N*4 (see Appendix A). 


is missing the ground radar депегаїес data (wind 


velocity), 


2. the radar wind speed and direction data appears but 


tne balloon data is missing, or 


р. all data is recorded (at the 7 or 8 mandatory levels 


оү 


Typically, either balloon or radar data is taken at 80 to 
100 points. 

Note that the Datain file and the program accept only 
conplete data sets. περ ανα пе ος 1999.97 for any data 
points that are missing. The program wili then do a semi- 
log interpolation for the missing values using one of two 


methods as follows; 
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2, 


1. if #RUNS = 1, all the missing values are assigned 
interpolated values. In this case, the program calcu- 
lates, plots and tabulates only one Cn? profile or 


2. if #RUNS = 2, the program rereads the raw datasets and 
interpolates only for missing values of WS and WD. Thus 
those datasets missing PRESS, TEMP and HUMID are deleted. 
In this case, a new Cn? profile is estimated, which is 
plotted and tabulated in addition to the previous one 
(1.6. when $RUNS = 1). Setting #RUNS = 2 nearly doubles 
the number of executions and thus the run time. The use 


of this option will be discussed in Chap. 4. 


The program may now be compiled (if necessary) and 
executed, it does not require batch processing and the 
output should be ready within several minutes. Note that it 
is possible to store many data sets in one's file, but the 
TURB2 EXEC reads from TURB DATAIN. Thus the user may find 
it necessary to either rename his/her DATAIN file or change 
the EXEC. 

If the user wishes to display a second Cn? plot simulta- 
neously with the model's estimated plot, the values for the 
second plot are stored in a second DATAIN file called 
REPEC? БАТАТ: This is not the same as setting #RUNS = 2. 
In this case a previously measured or estimated Cn? profile 
is input and the program merely plots the input values. А5 
in the previous case, the easiest method is to copy and 


rewrite an old file, however the format is as follows: 
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TABLE 3 
Reference Plot Datain File Format 


Card (s) Cols. Type Variable Remarks 





1 20 15 М number of data points 
to be plotted 


ΠῚ "eH Е10.4 Cn? reference values of 
Cne for Ponpe ΘΠ 
(not log values) 


" ο... AC T2022 hit absolute altitude 
μι to value 
О n 


А. COMPARISON 


Once the program was debugged and executing reliably, it 
was necessary to compare the output to that produced by the 
original program at NOAA. This procedure served as an 
acceptance test to validate the operation of the translated 
program. To this end, a random input and output dataset was 
reguested from NOAA. When the translated program was run 
with the test input on the NPS computer system, the output 
was found to agree closely with output produced by the orig- 
inal code executed on NOAA's CDC computer (see Tables 4 and 
5, Appendix D). Table 4 was produced using the first method 
described (i.e. #RUNS = 1) while Table 5 was produced using 
the second. 

As can be seen from Figure 4.1 and Table 4, when the 
first interpolation method was used (i.e. interpolated 
values are substituted for all missing datapoints), there 
are noticeable disagreements, but overall the I3M output is 
within 17.9% of the NOAA results (if one neglects the large 
Spikes). The difference in the accuracy between low and 
high altitudes was insignificant. There is one large excur- 
sion around 15 km. that is without explanation. Also, there 
were several Spikes in the NOAA data at 2 km. and at 6.3 km. 
that seem suspect and were not repeated in the IBM output. 

AS can be seen from Figure 4.2 and Table 5, when the 
second interpolation method was used (i.e. interpolated 
values are substituted for datasets missing WS and WD only, 
with the other datasets being deleted), the output tended to 
be in closer agreement with the NOAA results. However, in 


Several areas, the IBM data tended to respond more slowly 
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than the NOAA results (because of the fewer datasets) апа 
some of the spikes in the NOAA data were ‘smoothed' into 
insignificance. The underlying question is whether the 
resolution of the rawinsonde data is sufficient to predict 
the variations in the value of Cn? at the scale shown or 
conversely, is the data being overprocessed to the point 
that the spikes are merely 'noise'. To answer the question, 
it will be necessary to process more datasets and also to 
measure the Cn? profile experimentally so that comparisons 
can be made. 

Overall, with #RUNS = 2, the IBM output was found to be 
within 14.7% of the NOAA results. In the region below 12 
km., the agreement was within 20.1% of the reference data. 
The spike at 6575 meters seems to be produced by a glitch in 
the program (because the wind shear is zero for several 
consecutive slabs) and probably doesn't predict the actual 
Cn? accurately. Above 12 km., the agreement is excellent, 
averaging 5.0 and without any significant excursions. 

There are two primary reasons why the output doesn't 
agree exactly, as one might expect. The first is due to the 
fact that the data processing subroutine used in the NOAA 
program was not received with the rest of the code and tnus 
a substitute had to be written. The original (i.e.  NOAA) 
data processing routine was received later with the compar- 
ison test data. Although the two are very similar, there 
are some differences and thus slightly different values are 
input to the main program in each case. The second reason 
is that the NOAA program has undergone some modifications in 
the time between the receipt of the original code (fron 
which the IBM code derived) and the receipt of the compar- 
ison test data. While it is highly unlikely that any large 
differences would be introduced by the modifications, it is 
almost certain that some small dicferences crept in. Since 
the NOAA program itself is stili undergoing development, it 


is impractical to try to keep current translated copies. 
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However, with these differences in mind, the output of 
the translated code is quite satisfactory for the purposes 
specified and has the potential for being much superior to 
the accuracy of any of the previously available models. 
Also the ease of operation and useful graph/ tabular output 
format make it a very practical and useful addition to the 
NPS library. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

The NOAA program differs primarily in the smoothing 
algorithm used that produces a different dataset for the 
Main program. The translated program was modified so that a 
comparison could be made between the two methods of data 
processing described eariier. The intent of the comparison 
was to see if the first method was actually more accurate or 
if the larger number of datasets only tended to increase the 
'noise'. 

While the results indicated that the first method is 
better, it is difficult to make a judgement based on the 
results of a single dataset. However, for routine use, the 
first method is recommended (i.e. set #RUNS = 1 in the 
DATAIN file). 

One of the strengths of this model is that it requires 
only readily available rawinsonde data to be input, but that 
Simplicity is also a basic source of error that may occur in 
three forms. 

The first error comes about pecause the rawinsonde data, 
while available for a great many locations and times, is not 
available for all places and times. If one wishes to 
compare the model predictions with measurements from a radar 
Or Scintillometer, it is necessary to assume that the atmos- 
pheric parameters recorded by the rawinsonde are approxi- 
nateiy equivalent to those present at the measurement site. 
This may or may not be a reasonable assumption because of 
the spatial or temporal variation in the parameters. Гог 
Short distances (less 100 miles and without large variation 
in geographical £eatures like mountains or oceans), and 
times (2-4 hours depending on the time of day), the assump- 


tion is probably reasonable. Obviously, as distance or time 
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from the observation point increases, less correlation 
between the model and the measured data are expected. 

The second error derives from the fact that scintillom- 
eter measurements are time averaged line integrals, whereas 
the rawinsonde balloon carries an array of point sensors. 
Also the balloon rises over a one to two hour period along 
an obligue path as a function of the wind velocity, while 
the experimental measurements are usually made along a 
vertical one. For these last two reasons, even if the 
balloon launch site and the radar or scintillometer were 
co-located and readings were made simultaneously, there 
would still be some residual data errors. 

The third source of error occurs only if there are 
convective layers present in the atmosphere. Because of the 
difficulty in modeling the effects of such layers, the 
present model neglects then. 

With these limitations in mind, the model is still the 
best available and should be of utility to those whose 
research reguires the estimation of  Cn?(z). As the parent 
program continues to be improved, it may be worthwhile at 
some future time to add the corresponding updates to this 
version also. Additionally, i think that the program could 
be integrated into a propagation model, such ав G.U.T.S., 
for a significant increase in accuracy over the present 


model. 
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APPENDIX A 
TURB DATAIN FILE FORMAT EXAMPLE 


This appendix contains the dataset 
It is 


format required by the main 


used for the compar- 
an example of the 
Note that the file 


This is the format 


ison test. also to be used as 
progranm. 


contains two copies of the same dataset. 


3 7 


used when #RUNS = 2. If only one run is desired, the 
repeated dataset may de deleted (i.e. there would only be 
118 lines vice 236). 
ELOT? 1 (1-уе<, 2-по) 
#RUNS = 2 (1 or 2 only) 
# DATA SETS = 118 
PRESS. TSIP: HUMID. WS WD ALTITUDE 
(nbars) (GC) (5100) (n/s) (degrees) (neters) 
834.9 22.4 229.1 4.1 129; 16 10. 
828.0 2052 -261 DO 5) 1680. 
9 ЭОЕ 99979 Dal 155. 1930; 
778.0 15.5 5357 99959 29979 2200 
5.9 κ ο 95922 4.6 164. 2240. 
10518 0 12 Ὁ 2424 393959 7 2920 
9 275) 37 аа» 999.9 SE 172 2670. 
700.0 629 ΠΡ DEI 9 Ὁ 3100. 
Hu 9 5 999379 Oraal 176. 2140: 
5-9 οσο ο 90979 5 186. 55060. 
648.0 25 „863 0 5.9 3720. 
3999 999.9 SORS Deal 215. 9902 
615.0 2-21 2963 909.9 J997 4150. 
2993.9 mS G° 9 T 243. 4230. 
596.0 Hi 2 „690 99979 999.9 4390. 
Зоо J 3 Ον η 268. 460. 


1.0 
aod 
gon 
554.0 
озо 
52900 
SO 
513.0 
κ 
500.0 
399.9. 
883.0 
s 
460.0 
999219 
452.0 
58.0 
99959 
437.0 
93920 
410.0 
402.0 
400.0 
295959 
389.0 
ο 7 
99259 
92929 
349.0 
99529 
995.9 
316: 0 
95529 
300.0 
Qum 


- 714 
55 
99998 

- 664 
ооз 

2 2/2 
999739 

-411 
9998668 

se 
99979 

ου. 
99989 

. 054 
99929 

- 304 

„492 
ον 

- 204 
9999 

«207 

- 5955 

«17 
99865 

- 041 
9993 
29925 
59529 

. 242 
99979 
995.9 

.314 
ο 

„000 
399929 
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999.9 
E 
De 

999.9 
7.2 

999.9 
6.7 

999.9 
6.7 

999.9 
6.7 

999.9 
5.7 

999.9 
4.6 

999.9 

999.9 
4.6 

999.9 
5.1 
5.1 

999.9 

999.9 
4.1 

999.9 
3.6 
4.6 
4.6 

999.9 
5.7 
6.2 

999.9 
6.2 

999.9 
8.2 


99959 
279% 
278. 
99979 
274. 
99699 
276: 
99926 
280. 
9299 
295. 
96969 
266. 
999.9 
256. 
соо 
9E 
256. 
399959 
257: 
2102 
79953 
ЖӘНе 
278. 
99942 
266. 
214. 
270: 
9999 
259. 
264. 
ο 
211. 
9999 
280. 


4620. 
4680. 
4900. 
4960. 
5140. 
53209 
5410. 
5550. 
5690. 
5750. 
59804 
6010. 
6260. 
6380. 
6490. 
бә ТОЁ 
6640. 
6740. 
6760. 
6980. 
12305 
13708 
7410. 
75 10% 
7610. 
7800. 
8 120% 
8430. 
8560. 
8720. 
8990. 
9070: 
9230: 
9420. 
9470. 


DIIS 
99969 
283.0 
99929 
5-9 
399.9 
250.0 
299.9 
237.0 
Bag 
29979 
299.9 
98 Э 
9999 
200.0 
5.9 
lim». 0 
ου 
299.9 
29979 
29399 
150.0 
99 Send 
20929 
39979 
13020 
29529 
39949 
299.9 
117.0 
9 
108.0 
920969 
100.0 
39959 


-42.6 
5; 
= 19 
59929 
ЭЭ Э 
DE 
-46.0 
99539 
-48.5 
DS 5) 
99979 
99979 
ο. 
99979 
-5325 
53929 
2 
Зак 
999.9 
оа 
9929 
"σου 
299:9 
393929 
9222 
2596 
ο. 
ЭЭЭЭ 
999.9 
οσο Ὁ 
0 
-59.4 
999.9 
2595 
2535 9 


- 200 
99929 
- 000 
99929 
99929 
993929 
. 000 
29943 
-000 
99979 
99979 
gag 
99969 
999869 
„000 
595279 
„000 
999859 
299592 
999m9 
99979 
. 000 
99989 
OO oe 
99929 
. 000 
95259 
Заза 
DI 
- 000 
99929 
.000 
299.9 
„000 
DS) 
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99922 
‘Thiers 
3999 
14.4 
2125 
2322 
999533 
¿1.6 
995219 
24.7 
26.8 
Тез 
25819 
ҮЛЕН 
999.9 
D 
EX 
συ 
Leo 
2748 
239 2 
95523 
24.7 
24.2 
2m 
ον) 
24.7 
ош? 
225] 
99979 
20.6 
1920 
17. 0 
959.3 
16.0 


99979 
285% 
999:9 
283. 
ο. 
281. 
ο ο 
292% 
99979 
218; 
2 
2711. 
ο 
2з 
59929 
215. 
209. 
266. 
265. 
202. 
263. 
9° 929 
265. 
265. 
26 37 
59529 
262. 
264. 
264. 
39929 
264. 
2053 
259. 
999:9 
του. 


9510: 
2680. 
бо 
95 30 


10220. 
10500. 
10650. 
1077032 
11000. 
11020; 
11280. 
11540. 
117 9905 
120559 
12109. 
129202 
12600. 
12880. 
13160. 
13430. 
137 10. 
15930. 
13599051 
14280. 
14570. 
14830. 
14869. 
15160. 
15470. 
15500. 
15730. 
16000. 
16300. 
16480. 
16590. 


35 9 
999.9 
9929 
99979 
ολ 
DIES 
99979 
70.0 
ООО 
999.9 
999799 
62.0 
60.0 
3 
39969 
999:9 
99929 

50.0 
SE 
39959 
9395729 
42.0 
3293923 
39.0 
Зоо 

3060 
320399 
239229 
999379 

30.0 
99 9 
25650 


# DATA SETS 


9999 
gor 
999:9 
999.9 
999: 
39394,9 
Oo 5; 
-6 1.0 
go SS. 
3939 
9 
-58. 2 
295059 
999.9 
s 
999:9 
999739 
mee 
929.9 
339.9 
39959 
-54.3 
3299.9 
295220 
92929 
2928 
99979 
393953 
999.9 
=0 8. 2 
29079 
-47.0 
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99929 
999.2 
99999 
999.9 
99979 
999.9 
99979 
. 000 
999.9 
999.3 
99:98 

. 000 
.000 
gas 
9992 
39529 
93959 

. 000 
99985 
29) 
ЭЭЭЭ 

. 000 
goo 
«000 
95255 

. 000 
29923 
929929 
ο 

. 000 
99923 

- 000 


14.4 
13.9 
13.4 
12.4 
10.3 

9.8 
9.3 
999.9 
8.2 
ο 
5.7 
999.9 
4.6 
4.1 
4.6 
6.7 
4.6 
999.9 
4.1 
23 


-— 
° 
e 


О 
O 
Ne 
Ne 


О 
KO 
G © a 
ο 9 о ο 


a st à CO СУ ОО 
° е 
‘(кл л ІЛ іл ІЛ (ІЛ iD WM о СО 


268. 
264. 
269. 
270. 
266. 
266. 
266. 
299.9 
257: 
ΕΕ ΙΕ 
266. 
22122; 
262. 
266. 
290 
ӘРІГЕ 
31 z: 
999.9 
333; 
396. 
034. 
999.9 
064. 
ACT's 
347. 
99599 
007. 
238. 
229 
229: 
207. 
1355 


(delete here i£ #RUNS - 1) 
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16870; 
17150. 
17430. 
17 7407. 
1759505 
18260. 
18540. 
187 10. 
18820. 
19100. 
19380. 
19470. 
19670. 
19960. 
20240. 
20520. 
20500 
20830. 
21120. 
21440. 
21760. 
21950: 
220885 
22430. 
2277908 
22940. 
231105 
23450. 
23 T9 
20130. 
24430. 
24590. 


PRESS. 


(mbars) 


834.9 
828.0 
995,9 
778.0 
999.9 
750.0 
999.9 
700.0 
999.9 
999.9 
648.0 
999.9 
615.0 
999.9 
596.0 
999.9 
579.0 
999.9 
999.9 
554.0 
999.9 
529.0 
999.9 
513.0 
999.9 
500.0 
999.9 
482.0 
999.9 
460.0 
999.9 
452.0 
441.0 


TEMP. 
( C) 
22.4 
20.2 
992939 
ро 5 
99:519 
1250 
Τη 
ο 
99979 
9 =o 


HU MiD. 
(27100) 


«291 
- 28! 
ον 
ο. 
ago 
0424 
2999.9 
SOUS 
99929 
39949 
„883 
99959 
- 265% 
"πῃ 
«690 
299.9 
. 71% 
99559 
9999 
- 664 
95929 
220 
99979 
.411 
999379 
2619 
Θ᾽ 
ose 
Od 
. 058 
999.9 
. 304 
«192 
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WS 


(2/5) 


4.1 
эз О 
Oren 
99929 
4.6 
9 39579 
Diet 
399.2 
ӘЛЕН 
Do | 
29929 
6. ? 
99989 
Ta 
99529 
7.7 
99959 
του 
HE 
39959 
T 
299559 
Gay 
99959 
6.7 
99929 
6.7 
99979 
7 
993923 
4.6 
99959 
39922 


WD 


(degrees) 


120. 
99989 
158. 
93979 
164. 
399.9 
173; 
9 929 
176. 
186. 
39929 
212 
99959 
243. 
το ο) 
268. 
99929 
279. 
218. 
99959 
278. 
2939.9 
206. 
999.9 
2606 
99959 
2 15 
999;9 
266. 
995.9 
25б» 
99929 
9392. 


ALTITUDE 


(meters) 


1610. 
1680. 
1930. 
22098 
2240. 
25 20% 
20 70% 
3100. 
3140. 
351697 
3720: 
3930. 
4140. 
4230. 
4390. 
460. 
4620. 
680. 
4900. 
4960. 
5140. 
5 120: 
5410. 
9550. 
5690. 
5750: 
5980. 
6010. 
6260. 
6380. 
6490. 
6510: 
6640. 


9925 
437.0 
999.9 
410.0 
402.0 
400.0 
99965 
35920 
9059 
99959 
goo 9 
340.0 
99919 
399929 
316.0 
999.9 
300.0 
97989 
29020 
99. 
263550 
99959 
99979 
99925 
25020 
9959 
23 75 D 
99979 
999.9 
9559 
9999 
99329 
200.0 
99 ° 9 
112220 


99929 
- 205 
999.9 
-24.5 
προ ο 
22997 
99929 
EUM 
99979 
II 
ΠΗ 9, 
c 
99979 
99979 
209277 
99979 
-242.0 
999.9 
= 2.6 
99929 
-41.9 
99979 
599.9 
Jg gu 
-46.0 
ο 
21833 
5; 
9912533 
999.9 
99979 
899869 
25395 
ΠΡ 
-5979 


999.9 
«208 
η 
«287 
zs 
2417 
999735 
«081 
999.9 
999.9 
99923 
„242 
99999 
99929 
„314 
99329 
-000 
99939 
. 000 
99 29 
. 000 
2999 
99929 
QOEM 
2100010 
99929 
. 000 
9999 
999.9 
99959 
99979 
999229 
- 000 
99929 
- 000 
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4.6 
99948 
S 
921 
JOE 
9959 
4.1 
993253 
3.6 
4.6 
4.6 
JIII 
507 
6.2 
29959 
6.2 
99979 
9722 
99959 
11.5 
99979 
14.0 
2186 
23.2 
99909 
21.6 
99952 
2427 
26/5 
2ο 
2648 
2058 
399929 
2552 
ο 


256. 
999.9 
297. 
210. 
SS 
99979 
218. 
99359 
266. 
274. 
270. 
999 
259. 
264. 
η 
271. 
3999 
280. 
99949 
295- 
999.9 
ЕР 
zu 
ZEN 
99959 
282: 
98969 
219: 
LS 
Z7 Ve 
2 ПЕ 
273. 
99 
2735 
269: 


6740. 
6760. 
6980: 
7230. 
7370. 
7410. 
75 10" 
7616. 
7800. 
8 1208 
8430. 
85605 
8720. 
8990. 
900z 
9230% 
9420. 
9470. 
950: 
96808 
9810: 
9930. 
1022058 
10500. 
106502 
10770. 
11000. 
11020. 
112505 
11540. 
117985 
12050. 
12100. 
12320. 
12600. 


99939 
η 
999.9 
999.9 
1500 
9999 
5 
5-9 
130.0 
3 
399.9 
999.9 
117.0 
goo. 9 
108.0 
9.9 
190.0 
ου 
999.9 
2993.9 
299.9 
999.9 
299.9 
9993.9 
299.9 
70.0 
999.9 
299.9 
999.9 
62.0 
60.0 
Bus V9 
29979 
299.9 
9 


999. 9 
999. 9 
999.9 
999.9 
-56.3 
999.9 
999.9 
999.9 
-59.6 
999.9 
999.9 
999.9 
-58.0 
999.9 
-59.4 
999.9 
-58.3 
999.9 
999.9 
999.9 
999.9 
999.9 
999.9 
995.9 
999.9 
-61.0 
999.9 
999.9 
999.9 
-58.2 
-58.5 
999.9 
999.9 
999.9 
399.9 


9997 9 
93529 
SS >, 
999-9 

-000 
999.9 
999.9 
399.9 

„000 
99959 
999: 9 
99979 
. 000 
990929 

«000 
999.9 

. 000 
299.9 
999.9 
ЗЭЭ 
29979 
93929 
999213 
QS 
39923 
„000 
999.9 
ESOS 
99929 
. 000 

.000 
29979 
99973 
999,9 
999.9 


266. 
264. 
262. 
2686 
99979 
ӨЗ 
264. 
“τ 
99979 
262. 
264. 
264. 
93959 
264. 
265. 
2995 
993.9 
266. 
264. 
26%. 
DOE. 
ιν. 
266. 
266: 
266. 
53 
297. 
26 1. 
266. 
22929 
262. 
266. 
290. 
301. 
Su 


12880. 
13160. 
13430. 
137 10. 
13930: 
13990. 
14280. 
14570. 
14830. 
14860. 
151605 
15470. 
15500. 
15730. 
16000. 
16300. 
16480. 
16590. 
16870. 
17150. 
17430. 
177 10. 
17980. 
18260. 
18580. 
197103 
18820. 
Vo TOO. 
козо 
19470. 
19670. 
19960. 
20240. 
20520; 
20800. 


50.0 
ЭЭЭ 
99979 
999.9 

42.0 
323959 

33.0 
79929 

36.0 
999.9 
99959 
99969 

30.0 
999.9 

2:580 


=5 3:59 
ο ο 
999.9 
99929 
251. 
9995 
πο». Ὁ 
99929 
το. Ὁ 
III 
999 
99S 
=# 6 2 
39929 
- 7.0 


„000 
39979 
399.9 
99959 

- 000 
99979 

. 000 
ВЭЭ 

. 000 
999.9 
9932.9 
399929 

. 000 
99953 

. 000 
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`2 
Ne, 
“о 


— N ++ 
e п è 


“о 
ιο 
“о 
° 


uo 
uo 
woo 23 œ~ 
8 . 9 
{л л (л (л (n илло лоо ооа D 


еу =з „ ο ο O 
8 


99929 
333. 
356. 
034. 
ο 
064. 
267. 
347. 
99329 
O 0 
2395. 
229: 
229; 
207 
ιο» 


20830. 
2111286 
21440. 
21760: 
219508 
22080. 
22430. 
22710: 
22940. 
23110. 
23450. 
237 906 
24130. 
24430. 
24590. 


SUMMARY OF SUBROUTINES 


Subroutine name;  INTLGS 


Description; INTLGS integrates the functions ZL, FL, CNSQL, 
CDSQL, CQSQL and EPSLNL over the slab 


Inputs; L i C SHDDPIZUSHEAR, DODZ, Kl, CTHRSLB, ZLJ1IN, 
ZLMAX, ZLMEAN, SGMSA1, SGMST1 


το PL, PF; PCNSO, PCDSO, РСО5$0, PEPSLN 
Calls; SGMSHL, FNCTNL, GSSLGN 

Calied by; main progran 

Subroutine name; DATEKO 


Description; DATPRO interpolates for missing  rawinsonde 


data and calculates the weighted average and gradient values 
ШЫ; PRESS, TEMP, HUMID, 85, WD, 2, N, J, RUN 


Output; PARY, TARY, QARY, STBARY, SHEARY,  DQDZRY,  RIARY, 
THKARY, N 


Calls; INTER 
Called by; main progran 
Subroutine name; INTER 


Description; INTER performs a semi-log interpolation cf the 


input variables 
ИП ЕСЕ i B; C; D? NUMBER 
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Output; ANS 

Calls; none 

Called by;  DATPRO 
Subroutine name;  FNCINL 


Description; FNCTNL computes the value of the turbulence 
parameters (F, Cn?, and EPSILON) for a slab 


Inputs; P, T, Q, SIBLIY, SHEAR, DQDZ, RI, THKSES pK 
SGMSHR, SGMSTB, ZLMEAN 


Output; FL, CNSOL, СЕБЕ” СССР πρ τας 
Calis; BVRT 

Called by;  INTLGS 

Subroutine name;  BVRT 


Description; BVRT controls the computation of the turbu- 


lence parameters using a bivariate model 


inputs; Р, Т, О, >IBLTY; SHERR; ООШ τα ZL,  SGMSHPEE 
SGMST B 


Output; Е, CN5SO, CDS5SC COSO ЕРЕ РГ 
Cails; INTGSH 
Called рус иет 


{= ES one сш. Hun i ee “mamap. ЫЫЫ Е ӘӘ о. “sassa. “s ӘӘ ж Á s s q a COMES СӘНӘ НСАУ о s s “s жарар ӘР: ӘӘ - s a t r - “a “s < u w < “s 


Subroutine name; INTGSH 
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Description; INTGSH integrates the input functions over 


SHEAR using Gaussian quadrature 

Inputs; —SHRSLB, STBSLB, 21, SGHSHR, 5GMSIB 
output; F, CN2, BBS 

Calls; GSSLGN, FUNSHR, INTGS2, GSSLGR 
Called by;  BVRT 

Subroutine name;  FUNSHR 


Description; FUNSHR computes a function,  FTAU, from a 
Rice-Nakagami distribution of SHEAR's and a normal distribu- 
ШӘЛІ о STABILITY 


Inputs; TAU, TAUSLB 
Output; FRNTAU 
"eus; MBSS10 


Called by;  INTGSH 


Subroutine name; INTGS2 


Description; INTGS2 integrates the normal distribution 
function together with (SIABILITY)? times the normal distri- 


bution function over STABILITY using Gaussian quadrature 
Enputs: SIBRIC, STBSLB, SGMSTB 

ο σας CDFSTB, CN2STB 

Calls; STBINT, GSSLGR 


Called bv; INTGSH 


Subroutine name;  STBINT 


Description; STBINT integrates two input functions over 


STABILITY using Gaussian quadrature 

Inputs; STBMAI, STBHIN, STBSLB, SGMSTB, Nets, Sheik 
Output; CDFSTB;  CN2 TE 

Calls; GSSLGN, EXPNX 

Called by; INTSS2 

Subroutine name;  GSSLGN 


Description; GSSLGN calculates the abscissas and weights 


for Gaussian quadrature for limits from -1 to +1 
inputs; N 

Output; “ЛИ 

calls; none 

Called by; INTGSH, INTGLS 

Subroutine name;  GSSIGR 


Description: GSSLGR calculates the abscissas and weights 


for Gaussian quadrature for limits from zero to infinity 
Tnputs; N 

ο ασ τα 

Calis; none 


Called by; INTGSH, INTGS2 
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Subroutine; MBSS10 


Description; MBSS10 computes the value of the modified 
Bessel function of the first kind of order zero using the 
IMSL routine MMBSIO 


Inputs; X 

Output; ANSWER 

Calis; MBSSIO (IMSL routine) 
Called by;  FUNSHR 

Function name; EXPND 


Description; EXPND ccmputes the value of a normalized expo- 


nentiai distribution for ZL 
inputs;  ZL, ZLMEAN 

Output;  EXPND 

Calls; EXPNX 


Calied by; not in use presentiy, before modification, it 


was called in FNCTNL 
Function name; SHRSGM 


Description; SHRSGH computes the standard deviation of tne 


SHEAR distribution 

ποστς. -SGMCON, P, T, ALPHA, STBLTY 
τρις: 3SHRSGM 

Calls; none 


Cailed by; main program 
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Function name;  STBSGM 


Description; STBSGM computes the standard deviation of the 


fluctuating part of STABILITY 
Inputs;  SHRSGM, STBLIY 
Output; ΠΗ 

Calis; none 

Called by; main progran 
Function name; SGMSHL 


Description; SGMSHL computes the value of SGMSHL given 
SGMSH1 and ZL 


Inputs; TSGMSH1, ZL 
Output; SGHSHL 

Calls; none 

Called by; INTLGS 
Function name; PROB 19 


Description; PROB19 is intended to compute the probability 


function using polynomial approximation 
Inputs; “X 
Output; РЕОВ19 


Calls; none 
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Called by; not presently in use, before modification, it 
was called in INTGSH 


Function name;  EXPNX 


Description; EXPNX checks the input arguments to ensure 
they are within the operating limits of the computer. T 
they are not, it revalues them so that they are within 


limits to avoid underflows and overflows 
Inputs; X 

output; EXPNX 

Calls; none 


Failed by; STSINT, EXPFSH, EXPND 
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SLAB# 


OO ч OoN EU N = 


11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


A 
< 


28 
24 
25 
26 
2 
28 


1645.00 
1805.00 
2065.00 
2220.00 
2350200 
2595.00 
2885.00 
3129.00 
3350.00 
3640.00 
3825.00 
4035.00 
4185.00 
4310.00 
4425.00 
4540.00 
4650.00 
4790.00 
4930.00 
5050700 
5230.00 
5365.00 
5480.00 
5620.00 
5720.00 
9865200 
5995300 
6133.00 


МІРЛІТ (М) 
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SAMPLE TURB2 DATAOUT FILE (#RUNS = 1) 


CNe 


0.6429895076р-16 
0.8525603886р-16 
0. 1308713280р-16 
0. 1840405935D-16 
0. 1675688171р-16 
0-92089950-_ 95p= 17 
0. 9555537 126D- m 
0. 102851906 15-16 
0. 1845121879D- 17 
0.87456728887D-17 
0. 2248982453D-16 
0.70197468930-16 
0. 1676746736D-15 
0.207625 79 p= 15 
0- 1877315073D-15 
052996434 361D-16 
0. 28849 30028D-16 
0. 1875013105D- 16 
0.2 12356 103 Т Ше ЕБ 
0.25294 9526)- 16 
0.2375408888D-16 
0.1080201618D-16 
0. 1088445482D-16 
0.8812439910D- 17 
0- 9486586 7890517 
0.1078625192D-16 
0. 1777040099D-16 
0.2783679743D-16 


23 
30 
ВИ 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
БИ 
5/2 
5:3 
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56 
SI 
58 
59 
60 
61 
62 
63 


6320.00 
6435.00 
6500.00 
6575.00 
6690.00 
6750.00 
6870.00 
7105.00 
7300.00 
7390.00 
7460.00 
7560.00 
7705.00 
7969.00 
6° 75.00 
8495.00 
8640.00 
8855.00 
9030.00 
9150.00 
9325.00 
9445.00 
9490.00 
259 3:00 
9745.00 
9870.00 
10075.00 
10360.00 
10575.00 
10710.00 
10885.00 
11010.00 
11150.00 
11410.00 
11665.00 


o 


0253856 149385D- 16 
0.2373786 114D-16 
0.2249937021D- 16 
0.2985857815D-18 
0.1865800713D-16 
0.1937612712D-16 
0.1052763186D-16 
0.1490718353D-16 
0. 1813311719D- 16 
0.28923847300-16 
0.2745028922D- 16 
0.1685917667D-16 
0. 1646131759D- 16 
0. 163875652523D-16 
0-5595951535D=17 
0-1192759561Dp-16 
0. 1121599474D-16 
0:96 7/16908 70D-17 
0.10450687732-16 
0:12528298310-16 
0.3279598721D-16 
0.3337614124D-16 
0. 1042557235D- 15 
0.65044801240-16 
0.4285484888D-16 
9-8093830537D=16 
O23 11097 53455D= 15 
0. 7411441978D-16 
02,027 2454 522D-16 
0. 1333465354D-16 
0.583 1460 262D- 16 
0.5737078697D-16 
0. 3580049 158р-1є 
0. 1080 384642р-16 
0295339098 12:7D- 1 7 


6 ἡ 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
{ 
76 
77 
7δ8 
19 
80 
81 
32 
83 
84 
85 
86 
87 
88 
89 
90 
971 
22 
93 
94 


26 
mT 
98 


11920-00 
12075.00 
12210.00 
12460.00 
12740.00 
13020.00 
13295.00 
13570.00 
13820.00 
13960.00 
14135.00 
14525.00 
14700.00 
14 845.00 
15010.00 
15315.00 
15485.00 
15615.00 
15865.00 
16 15 0.00 
16390.00 
16535.00 
167 30.00 
17010.00 
17290.00 
17570.00 
1784 5.00 
18120.00 
18400.00 
18625.00 
1876 5.00 
18960.00 
19249.00 
19425.00 
19570.00 
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0. 1521933891D-16 
0. 1874242331D- 16 
0.1619829771D-16 
0.23744955062-16 
0.25981012840-16 
0.15491311822-16 
0. 1054 737223D-16 
0. 1540170661D-16 
0.1251849086D-16 
0. 1213451593р- 16 
0.1147173704D-16 
0.5974197139D-18 
0.1379922626D-17 
0. 1774158125D-17 
0. 28833704 86D-17 
0. 1067950577D-16 
0.1216896389D-16 
0.9765411422D-17 
0.8095732044D-17 
0. 1096763098D-16 
0.9651295968D-17 
0. 9 2076335030- 17 
0.8188833993D-17 
0.6958387906D-17 
0.7090243292р-17 
0. 6583973759D-17 
0. 7582720669D-17 
0.5836168579D-17 
0.4927868104D-17 
0. 5680862268D-17 
0.5600901504р-17 
0. 4850331731р-17 
0.7716729501D-17 
0.6594158104D-17 
0.4023897304D-17 


99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
uo 
110 
111 
112 
113 
114 
115 
116 
117 


551 51018 
20100.00 
20380.00 
2066 0.00 
20815.00 
20975500 
21280.00 
21600.00 
21855.00 
22015.00 
22255.00 
22600.00 
2285500 
23025.00 
23280.00 
23620.00 
23900200 
24280.00 
24510.00 
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0.3734642299р- 17 
0.4261320896р-17 
0. 4530608595р-17 
0.5746324061D-17 
0.4799269686D- 17 
0.4589887089D-17 
0. 06832589 6020-17 
0.2859715419D-17 
012551690 196D-17 
0.3095974497D-17 
0. 3229411536D-17 
0. 2 199815305D-17 
0.2342476571D-17 
03 2279193328D-17 
0. 2015465644р- 17 
0.24060930328D-17 
0-2365552736D-17 
0.1883079643D-17 
0. 1813421204D- 17 


Table 4 
obtained from 


Specified altitudes. 


CCBPAR 


the NOAA 


mman m. Qum qam ms 


ISON OF 


program and 


RESULTS 


illustrates the corresponding values 
the IBM 


The fourth column lists 


nolel at 


of Cn? (z) 
the 


the percent 


difference in the IBM value relative to the NOAA value. “he 
#RUNS was set at 1 which resulted in the MIDALT values being 
the same for both. 


MIDAIT 


1810.00 
2070700 
ο 
238000 
2690.00 
2880.00 
3120200 
335 0:00 
3640.09 
3830.00 
4040.00 
4190.00 
4319.00 
4420.00 
4549.00 
4650.00 


(m) 


NOAA Cn? 


0:7906F 16 
0.1029Е-16 
951555 20 
О-5120Е—20 
0.108 1Е- 16 
0291278817 
0.140SE- 18 
O0.2106F=-15 
0.37 5 p= 15 
0.2510EF-16 
0555 ΟΞ“ 
О„1529Е=15 
0.2080F=15 
02207345515 
0-28 210216 
0.3446E-16 


TABLE 4 
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Comparison of Program Output with $RUNS 


MODEL Cn? 


ο SEED 
0. 1309&- T6 
0- 1880EËE-16 
0. 1676E-16 
0.929051 
0- 9556E-17 
0. T022 EF T6 
0. 1845Е-17 
ο, δη 
0. 2249Ε-16 
0ο. 9205-16 
0. 1677£- 15 
OSZ20768— 15 
0. 1877E-T5 
0. 29965711E 
OZ EIE 


7.812 
27.185 
354574 .06% 
327182 .8 1h 
14.06% 
4.707 
7199.647 
776.132 
220620105 
10.40% 
24.02% 
9.66% 
0.18% 
ΑΞ; 
4.04% 
16.282 


790.00 
930.00 
5050.00 
9230.00 
5360.00 
5480.00 
5620.00 
5720.00 
5870.00 
€000.00 
6140.00 
6320.00 
6430.00 
6500.00 
6570. 00 
6690.00 
€730.00 
6870.00 
7190.00 
7300.00 
7390.00 
7460.00 
7560.00 
7700.00 
7960.00 
8280.00 
8500.00 
8640.00 
8850.90 
9030.00 
9150.00 
9330.00 
9440.00 
9499.00 
9600.00 


Oa 26816 
0.194 2E- 16 
0.2852E-16 
ОФ 2 2885516 
0-1090E-16 
9 ΠΕΙ 
0 78/98 37-17 
079317 Στ 
Оос EO 
Qo T9r-TO 
02270 Sere 
029 IESO 
0.3885E-20 
0.3467E-20 
0. 1165F8=20 
0: 1912Е-16 
О- 188 7Е- 16 
0. 106 92-16 
0.146 2E-16 
ОЕЕО 
0.1974Е- 16 
ο πο ЕЕС 
ο σοι »- 16 
US SOON NM. 
ОЭ Ти 
0.6828Е-17 
02105355F- T6 
029070 15-17 
Qo 9235—17 
05392 Sr 9 
О- КУЗ 5Е-— 16 
2-29 2Е 16 
0.3479E-16 
0.9 T7 57-6 
sg 2E- 6 


©] 


OQ RETO SES T6 
Q. 2126E=16 
0. 2529Е-16 
OB 2331 Sits 
ПЕ ШЕШЕ = 16 
0. 1088F=16 
02 88 12E=17 
0.9487Е-17 
0; 1079J9E-16 
0. 1777E-16 
0.2784E-16 
0. 3 05/5 8:16 
0.2374Е-16 
0.2250E- 16 
022996ES168 
0. 1866E-16 
0. 1938E-16 
0: 1055516 
0.1491Е-16 
0. TONS E- 16 
02128902516 
0.278%85E-16 
0.1686Е-16 
0. 1640E-16 
0. 1689/55116 
0669 p= g 
ο ilo 32 STO 
ΙΙ 518 
ο. οσο) 
0. 1085E- 16 
ΤΕΙ: 
0. 3280E- 16 
0. 33939516 
0. 10438-15 
0.6504Е-16 


8.637 
9.45% 
ο. τση 
5.62 
0.90% 
5.78% 
1.90% 
ο ο x 
1.475 
9.76% 
> 
σος 
6810483. 3 Ίο 
6538857 «373 
25463.847 
2.97% 
2.685 
oes 
1 
5.0 35 
46.52% 
51.49% 
4.39% 
dM 
82.85% 
Ὁ 
13:06% 
28.90% 
66.095 
76.265 
1.86% 
7.85% 
4.067 
T3763» 
11.91% 


9750.00 
9870.00 
10070.00 
10360.00 
10570.00 
10710.00 
10880.00 
11010.00 
11150.00 
11410.00 
11670.00 
11920.00 
12080.00 
12210.00 
12460.00 
12740.00 
13020.00 
13290.00 
13570.00 
13820.00 
13960.00 
14140.00 
14430.00 
14700.00 
14850.00 
15010.00 
15310.00 
15480.00 
15610.00 
15870.00 
16150.00 
16390.00 
16540. 00 
16730.00 


0.4248 0E- 16 
0.535015 
0.528 4E- 15 
0. T69 9 ES 
0.16255-16 
ОТЕ Ee 
0.554 6E- 16 
ο οἱ. RENI 
0.357256 
0. 125 8 E205 
0.11135= 16 
0.1192E=16 
O. 1600EF- 16 
0.1566Е- 16 
0.2469Е- 16 
0.2574Е-16 
ο το της 
ОО ЕИБ 
0.138 7Е- 16 
0.1048Е- 16 
0.8874E-17 
(129222517 
0292855 07 
0.500 35-17 
DIS τ; 
0.1034 E- 16 
029622537 
0. 1059 ESI 
ο ο ορ 
0. TOS Su 
0 DOSSE-N6 
ο σετ. 
ο. 
0-7117F-17 
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0.4285E-16 


0. {0555516 
0.111.515 
0. 1&8 11E- 16 
0. 12425515 
0. 139 3BE- 16 
0.-5831E-16 
0.573765 16 
0.353015 16 
0. 1080E- 16 
02 953^ Ed 
ο. 1522Е=16 
0. 1874 E- 16 
0. 1620E- 16 
0.23 25 E26 
0. 25985 = 16 
0. 1549Е- 16 
0. 1055E-16 
0. 1580E-16 
0 20025 252316 
0. 1213E-16 
О ЕЕЕ 16 
0-5978E=-18 
0. 1380E-17 
0. 17748 E- 17 
0. 2953] 17 
0. 1068E- 16 
ο. ΕΖ Ε 16 
0-25-17 
0.8096Е-17 
0. 109 7E- 16 
0. 9639527 
0.9208E-17 
0. 8 189E- W7 


1.07% 
23. 38 
41. 135 
16. 93% 
21.70% 
1.71% 
5. 1 
10.503 
0. 23 % 
16. 12% 
14. 34% 
27. 68% 
17. 19 5 
3.44% 
3. Sam 
0.94% 
3.12% 
8. 80% 
11.04% 
19. 455 
36. 74 % 
39:529 
932515 
82.76% 
81.857 
72.115 
23. 86% 
14.915 
5„ 20 
3. 945 
3. 18 
5.51 
10. 74% 
15.06% 


17010.00 
17290.00 
17570.00 
17840.00 
18120.00 
18400.00 
18620.00 
18760.00 
18960.00 
19240.00 
19420.00 
19570.00 
19810.00 
20 100.00 
20380.00 
20660.00 
20820.00 
20970.00 
21280.00 
21600.00 
21850.00 
22020.00 
50-00 
22600.09 
22860.00 
2:059. 00 
23280.00 
23620.00 
23960.00 
24280.00 


ο δε 
0 0150p-17 
025069 27-207 
U26798E-17 
0.4908E- 17 
0.466 1Е- 17 
ums 10 E= 
0m»5967E- 17 
ο ο πας 
0.640 1Е- 17 
0.4795Е-17 
95 ЗОЗӘРЕІ ТІ 
mr 362r-17 
0.B735E-17 
0.4833E- 17 
0538055 17 
033779E-17 
Cis 1226817 
ο ολες] 
ΑΡΑ 8, 
ο το RET 
C202 te 17 
02:3: 1072: E3057 
ο Io Or 
ο ο Επ 
Dp SE 7 
02235 CEI? 
gno DOSE ТІ? 
DS WEST 
0. 1840E-17 
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09609568E- 17 
ШЕ 7090Е- 17 
0. 6584E- 17 
0. 75893Bp=17 
0" 556306E- 17 
ο θες 
0a 569E- 17 
0. 5601E- 17 
0. 4850 E- 17 
ОЛЕ ТИЕУ 
0.6594Е- 17 
0.4024E-17 
Ges 7S5E-17 
O.4261E-17 
0. 4531E-17 
0.5746Е-17 
0- 4799E-17 
ο ο ΟΕ 
0. 4683E-17 
De b) 
255292017 
D 30965E-— 17 
νυ 9p= 17 
0. 2200E- 17 
0. 2342E- 17 
ОСОБЕ 
υ 5155; 
0.2861Е-17 
ο ο. 


0- 18832-17 


15.204 
19-225 
15.87% 
11.54% 
18.91% 
Serr ai 
9.04% 
4.547 
3. 99% 
20.56% 
37525 
4. 84% 
14. 385 
10. 00% 
6. 35% 
Ооз 
27:00. 
47.02% 
36. 86% 
1. 26% 
0. 78% 
2.48% 
1. 81% 
2599 
Veena 
5.629 
14. 53% 
13.467 
13: 65% 
2.34% 


Table 5 illustrates the corresponding values of Cn? (z) 
obtained from the NOAA program and the IBM model at the 
specified altitudes. The fifth column lists the percent 
difference between the IBM value and the interpolated (seni- 
log) value from the NOAA data. Note that the different 
interpolation method (i.e. #RUNS = 2) produced different 
values for the MIDALT's thus making direct comparison 


inpossible. 


TABLE 5 
Comparison of Program Output with #RUNS = 2 


NOAA data IBM data % Difference 
Cn? Alt. (m) Cn? Alt. (m) 

26430E- 16 1645. 
.7908Е- 16 1810. 

429m 16 1940. 47.8 
-1029Е-16 2070. 

.1671Е-16 2360. 5727 
.1081E-16 2600. 

92 5p 17 2810. 208 
.91272-17 2880. 
. 1409Е- 18 3120: 
„2106Е-18 33508 

.2361E- 18 3410. 4.2 
ο 5 18 3640. 
20758195 — 16 5550: 

209021516 оло 7 
. 5660E- 16 040. 
-15298-15 4190. 

.1777Е-15 4265. ο 3 
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ESXUDSÜE- 15 
870 7 YE- 15 


ESO 'E- T6 
~-3446E-16 
ИР 2СЕ- 16 
SSH 2rZ- T6 
ESI 52E-16 


2 ӘВЕ- 16 
EU SÜE- 16 


ШУ? JE-16 
298 3E- 17 


E17 7E- 17 
ШОЗЕ- 16 


1155-16 
a270 5E- 16 


E 39E- 16 


ШІБІ2Е-16 


2155 75- 16 
mure T6 


ІН 6 2Е- 16 
BI 2 0E- T6 
EIS T HE- 16 
. 15 12E=S16 


6 516 


4310. 
4420. 


540. 
650. 
4790. 
4930. 
5050. 


5230. 
3360. 


5480. 
5620. 


5720. 
5870. 


ου 


6150. 


63 20; 


6590. 


6750. 
6870. 


7100. 
7300. 
7390. 
7%60. 


7300: 


165-515 


. 1777 £- 16 


ο ο. μας 


. 1068E- 16 


.8406E- 17 


ΟΤΕ 6 


22070 b= 16 


‚2362Е=16 


ЛЕ 18 


EU e END 


.1142z- 16 


«1920 E- 16 
52057-16 


z 1563P 16 
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4505. 


4790. 


5140. 


5435. 


56505 


5550. 


6155. 


6445. 


6575. 


67007 


69953 


7300. 
73907 


7510. 


S 3 


7.8 
9526 


11.3 


ο 
46.3 


29099 E 
292 Em 


„6828Е=17] 
571402510» (5 
ο.” 


οτι 
ο ο 
«1555-15 


«310 9E dO 
. 3479E- 16 


429128970 


- 4240E-16 
223:3/5 015215 


; 169 9916 
СЕБЕ το 
s13 ES Бо 


.5516E-16 
EM asc G 
. 390 25b 
41209-16 


«111515 :16 
Pei 116 
«16005- 16 
.1566Е-16 


:215069EFE=16 
+29/4E5 16 


7700. 
7960. 


8280. 
8500. 
8640. 


8850. 
9030. 
9150. 


35 30. 
9440. 


9600. 


9750. 


10360. 
10570. 
10710. 


10580. 
11010. 
11150. 
11410. 


11670. 
1-495720; 
12080. 
12210; 


12460. 
12740. 


ιδεα 


οσο 


“οσους 16 


23837 5- 16 


2:982 5B- 16 


.1586E- 16 
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